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Abstract 

We study eta meson properties in the infinite nuclear matter and in atomic nuclei 
with an emphasis on effects of the eta coupling to iV*(1535)-nucleon-hole modes. 
The A r *(1535) resonance, which dominates the low-energy eta-nucleon scattering, 
can be seen as a chiral partner of the nucleon. The change of the chiral mass gap 
between the N* and the nucleon in a nuclear medium has an impact on the proper- 
ties of the eta-nucleus system. If the iV*-nucleon mass gap decreases with a density 
increase (chiral symmetry restoration) the calculations show the existence of the 
resonance state at the energy about 60 MeV and two bound eta-nucleus states 
with the binding energies about —80 MeV. These states can have strong effect on 
predicted cross sections of the 12 C(7,p) 11 B reaction with eta-meson production. 
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1 Introduction 

The eta bound states in nuclei were theoretically predicted by Haider and 
Liu [1J. This prediction relied on the first estimation of the rjN interaction 
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just made that time by Bhalerao and Liu [2|, who obtained an attractive r]N 
scattering length. Being electrically neutral the eta meson can be bound only 
by strong interactions. Simultaneously, the strong interactions induce nuclear 
absorption of the bound eta meson, which wave function has to have large over- 
lap with a nucleus. As the result the position and the width of the bound level 
are determined by a complicated interplay of the nuclear density distribution 
and the strength of the in-medium r)N interaction. Therefore the experimen- 
tal proof of the existence of eta bound states could provide information on 
modification of meson properties in nuclear matter complimentary to that ob- 
tained from pionic and kaonic atoms [3J. Particularly, the existence of such 
states in light nuclei would set important constraints [4]. Currently, various 
theoretical calculations predict the natural widths of eta bound states to be 
larger than the level spacing [5|6|7f8|9|10flllll2j . This could explain problems 
in experimental identification of the eta bound states. 

The first experimental search for eta bound states in (tc + , p) reactions [13] gave 
a negative result as being aimed only at the narrow states as it was proposed 
in Ref. pp. At the same time the measured excitation functions for reactions 
d(p, 3 B.e)r] [14] and 18 (tt + , tt~) 18 Ne [32] showed some enhancement near 
the r) production threshold which could be a hint for a bound state. Sokol 
and Triasuchev argued in [16] that the eta bound state can be more easily 
identified if one tags on the iV(1535) resonance in the final state by observing 
its 7T + n decay. Such an identification method was realized in experiments at 
the Lebedev Physical Institute studying (7,p) reaction on the 12 C target. 
The observed enhancement of n + n production occurring when the photon 
energy exceeds the eta production threshold can be interpreted as a signal 
of //-nucleus formation [17J. Observation of an eta bound state in 3 He with 
the binding energy (-4.4 ± 4.2) MeV and the full width (25.6 ± 6.1) MeV in 
photoproduction reactions at MAMI was reported in [TS|. The bound state 
was identified by registration of two decay modes with an explicit eta in the 
final state and with n°p. Interpretation of these experimental results in terms 
of the eta bound state was, however, questioned in [19]. A new support for 
?7 3 He bound state arises from the recent study of dp — > 3 He 7/ reaction at 
COSY [20]. Other experiments have been proposed for searching eta bound 

Proton pickup processes are powerful experimental tools for studying light 
mesons in nuclei. They allow for the recoilless creation of a meson inside a 
nucleus by fine tuning the energy of an incident particle. It is well-known, 
e.g., that the (d, 3 He) reactions on nuclear targets lead to formation and iden- 
tification of the deeply bound pionic atoms [2"4"|25f26|l2"T] . The (7,2?) and (n,p) 
reactions have also been proposed as good experimental tools for formation 
of eta-nucleus systems. In the recoilless production, the eta meson is created 
almost at rest and may be absorbed within the nucleus. Although the eta 
meson is not directly observed in the final states, we can extract properties 
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of the eta-nucleus system by measuring of the production cross section as a 
functions of the momentum of the emitted particle [T1I10|I11|I12| . Thereby one 
has access to the energies corresponding to eta mesons bound in nuclei. This 
is one of the differences from the eta production reactions on nuclear targets, 
where the eta meson emitted from the nucleus is observed [28ll29f30| . 

The prominent feature of the eta meson is that the eta-nucleon system cou- 
ples dominantly to the iV*(1535) (Su) resonance at threshold energies [31]. 
The r]NN coupling is, on the contrary, much smaller than rjNN* and irNN 
ones [32J and is basically irrelevant for analyses of the NN scattering and the 
threshold 77 production on a nucleon [21I33] . Therefore the low- momentum eta 
meson would dominantly excite iV*-nucleon-hole (iV* — h) states in nuclear 
medium. The difference between the eta mass and the N* — N mass gap is just 
50 MeV, which is smaller than the N* width T^* ~ 75 MeV and is comparable 
with the depth of the attractive potential acting on the nucleon in a nucleus, 
~ 50 — 60 MeV. Hence, one could expect the eta meson nuclear dynamics to 
be sensitive to modifications of nucleon and N* properties in medium. 

The question about in-medium dynamics of the iV(1535) resonance addresses 
an interesting aspects of the chiral symmetry of strong interactions. Being 
the lowest lying baryon with the opposite parity to the nucleon, the iV(1535) 
can be viewed as its chiral partner [34,35,36,37,38]. Analogously one considers 
pairs (cr, 7r) and (p, ai). If now the chiral symmetry is getting restored with 
an increasing baryon density the parity partners will become degenerate. This 
picture has several interesting phenomenological consequences [3"4"|35f38| . The 
important implication for the eta meson physics is that the N — N* mass gap 
decreases in nuclear medium and at some density goes below the rj mass [391136] . 
This changes the sign of the rj meson optical potential making it repulsive at 
saturation density as found in [rU|ll||12j . 

The parity-doublet concept corresponds to a linear realization of chiral sym- 
metry in terms of hadrons. In a non-linear realization the parity partner is 
not present in the Lagrangian as an independent degree of freedom but rather 
manifests itself as a pole of the scattering amplitude generated dynamically 
by coupled channels. So, the iV*(1535) can be viewed as a KH and KA quasi- 
bound state [4"0~|41ll42|r4"3"] . In general the dynamically generated states could 
be very sensitive to in medium modifications of the meson-baryon loops they 
are made of. For instance the dynamically generated A (1405) resonance would 
dissolve or become very broad in nuclear medium [4*4l45ll46|T4T] . In the case of 
the iV*(1535) the in-medium modification of kaon-hyperon loops is expected to 
be smaller [4"4"f4"8"] since there is no Pauli blocking for hyperons and K mesons 
interact rather weakly with nucleons. However, even in this case both hyper- 
ons and nucleons can couple to an attractive scalar nuclear mean field, which 
is responsible for a decrease of baryon masses in nuclear matter. Thus, the 
shift of the N* with respect to the nucleon depends on how strongly hyperons 
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and nucleons couple to the scalar field. If the couplings are equal then the N 
and N* masses decrease roughly in the same manner that keeps the N* — N 
mass gap constant. It is known, however, from physics of hypernuclei that the 
hyperons are bound weaker than nucleons. This could lead to an increase the 
N* — N mass gap in the nuclear medium. 

As discussed above, these two different pictures for iV(1535) imply distinct 
phenomenological consequences for eta meson dynamics in the nuclear mat- 
ter [TTp2] . In the parity-doublet model the eta optical potential is shallower, 
that reduces the binding [T01I11II12] . 




Typically, calculations done in previous works operated with energy-indepen- 
dent optical potentials. In the case when iV* — iV gap becomes close to or 
smaller than the eta mass it is not a reliable approximation anymore, and 
the energy dependence must be taken into account explicitly. The eta meson 
quantum numbers can be carried now not only by the eta-particle modes 
but also by the iV*-nucleon-hole modes. In the present paper we study the 
propagation of these modes in nuclear matter and finite nuclei. We discuss a 
possibility of the formation of an "eta nucleus" by iV* — h modeJ"*"! 

In Section 2 we construct the eta meson Green function in nuclear matter ex- 
ploiting the iV* dominance in the r]N scattering. We calculate the spectrum of 
eta excitations in the infinite nuclear system and analyze how the eta spectral 
density changes with the nuclear matter density. Section 3 is devoted to eta 
excitations in finite nuclei. First we discuss the spectral function convoluted 
with the density distribution in n B. This gives an approximate view on how 
the nucleus would respond to the external source with the eta quantum num- 
bers. For better understanding the observed bump structures we solved the 
Lippmann-Schwinger equation for the eta-nucleus scattering amplitude in the 
momentum representation and identify poles of the amplitude on the com- 
plex energy plane. In Section 4 we calculate missing mass spectra in the (7,p) 
reaction on 12 C target associated with the eta meson formation. 



2 Eta meson in nuclear matter 

2.1 Eta self-energy in nuclear matter 

In this section we investigate the eta meson in infinite uniform nuclear medium. 
The propagation of the eta meson is described in terms of the in-medium Green 



Similar interpretation of the "eta nucleus" in terms of both eta-particle and 
resonance- hole modes was proposed in [IB] . 
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Fig. 1. Diagrams for the r/ self-energy in the nuclear matter, (a) and (b) are the 
direct and crossing diagrams, respectively. 

function with the self-energy H v (u, k] p) 

1 

G v (u,k;p) 



uj 2 — k 2 — m 2 — n^(o;, k; p) + ie 

Here uj and k denote the eta energy and momentum, p is the nuclear density, 
and m v is the vacuum eta mass. We consider the isospin symmetrical nuclear 
matter. 

The eta nucleon interaction is dominated by the N* resonance which con- 
tributes to the s-wave scattering amplitude. The higher partial waves are less 
important for slow r) mesons created in reactions at recoilless conditions. Thus, 
the main contribution to the eta self-energy is due to the excitations of N*- 
nucleon-hole modes as shown in Fig. UJ Exactly these diagrams induce the 
strong energy dependence of the self-energy. 

The T]NN* coupling is defined by the Lagrangian 



C vNN *=g(NN* + N*N)r ] . 



where the coupling g can be estimated as g ~ 2 from the vacuum N* — > rjN 
decay width 



2 ("^* + m N? ~ m l f, 2 2 W 2 T~2 2 

T N *^ vN = g — 3 L J{m z N , - m z N - m 2 ) 2 - Am%m 2 

10 71-77%* v ' ' 



with the vacuum values, Tn*-*^n — 75 MeV, tjin = 939 MeV, m^v* 
1535 MeV and m v = 547 MeV. 

In the gas approximation the self-energy depicted in Fig. [T] reads 



9 2 P 



H„(uj, k: p) = 

" w- (mUP )-m* N {p ))-^ + iT N ,(u J ;p) 2 
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where the first and second terms in Eq. ([3]) correspond to the diagrams (a) 
and (b) in Fig. [TJ respectively. We neglect here momenta of the nucleon and 
the N*. The linear density dependence in the numerator of Eq. ([3]) is obtained 
by the Fermi momentum integration over the nucleon hole state. The density- 
dependent in-medium masses of the N* and iV are denoted (p) and 
m* N (p), respectively, r^. (u; p) is the in-medium N* width determined by the 
iV* self-energy, T N *(u; p) = —2 QSjv«(m^(p) +u;p) . There is no width in the 
second term in Eq. ([3]) since QT lN *(m* N (p) — cu; p) = 0. 

Below we use the chiral doublet model [38] for the medium modification of 
the TV and N* masses. In this model, the mass difference of the nucleon and 
iV* is determined by the chiral condensate (a) as 

m* N .(p)-m* N (p) = A((T) (4) 

with the parameter A = {tun — tun*)/ (v)o, where (<r)o is the vacuum chiral 
condensate. Assuming a partial restoration of chiral symmetry in the nuclear 
medium and parameterizing of the sigma condensate as a function of the 
nuclear density [49J, (a) = (1 — Cp/p )(a)o, we obtain the in-medium mass 
difference 

m* N *{p) - m* N (p) = (m N * - m N )(l - C p/p ) . (5) 

The parameter C represents the strength of chiral restoration at the normal 
density and its reasonable values are in the range of C — 0.1-0.3. The value of 
the C parameter can be evaluated within the chiral doublet model providing 
C ~ 0.22, which is read off the density dependence of the chiral condensate 
calculated in Ref. [39]. In this paper we use C = 0.2 [TU|ll||12j . The mass of 
the eta meson is assumed to be constant in the nuclear medium [5T)llo"Tj (no 
scalar mean-field potential). 

There are several contributions to the imaginary part of the iV* self-energy, 
among which the main modes are the pionic decay, N* — > ttN, and the two- 
body absorption, iV*iV — > ttNN [6]TT0] . The decay process iV* — > i]N is for- 
bidden in nuclear matter due to the Pauli blocking. The width of N* — > nN 
reaction calculated in Ref. |39j within an improved parity-doubled model is 
depicted in Fig. [2] (dashed lines) as a function of the energy for different nu- 
cleon densities. The steep rise of the width is due to the pseudoscalar ttN*N 
coupling motivated by the chiral symmetry. For C = 0.2 this partial width is 
decreasing with a growing density as the chiral symmetry is getting restored 
and pions decouple from baryons. For C = 0.0 the width T n *^ tN does not 
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Fig. 2. Energy dependence of the N* widths: the partial width Tn*^ttN (dashed 
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a v N [fm] 
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a v N [fm] 
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0.55 ± 0.02 + i 0.3 


0.88 + i 0.27 


0.98 + 2 0.37 


Ref. 
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[53] 


m 


[55] 



Table 1 

Summary of the values of the r]N scattering length. 

depend on the density. The N*N — > tcNN reaction has been calculated in 
Ref. (TO]. The total N* width is depicted in Fig. [2] by solid lines for C = 0.2 
and by dash-dot-dot lines for C = 0.0. The analysis done in this section for 
the case of the homogeneous nuclear matter shows that, in principle, it is suf- 
ficient to take density and energy independent value for the total width T N * 
given by the vacuum value TN*->nN = 75 MeV. Later in the calculation of 
the formation spectra of the eta-nucleus system, we will use more realistic N* 
decay widths reported in Refs. [T0|llf58j . 

We should verify how our model describes the low energy r)N interaction. The 
r/N scattering amplitude can be related to the r\ self-energy as T vN (u,k) = 



— dH v (u, k; p)/dp 
length 



p=0 



Using this relation, we calculate the r/N scattering 



= T vN (m v ,0) ^ Q 55 . Q 3g fa 
v Arr(l + m v /m N ) v ; 

Note that T vN (u) = m^) = 0, since the N* — > r/N channel is closed at the 
threshold. The values reported in the literature are summarized in Tabled] Our 
value of the scattering length is consistent with the other model calculations. 
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Fig. 3. The real (left panel) and imaginary (middle panel) parts of the solutions, 
u>i, of eq. ([8]) and the occupations factors ([9]) (right panel) for vanishing momentum 
as functions of the baryon density. The dashed lines on the left panel denote the 
eta and N*-h modes without the channel mixing (g = 0). Calculations are done for 
C = 0.2. 

With the positive rjN scattering length, the nN interaction is attractive at 
the threshold. Thus, the self-energy of the eta meson is also attractive at low 
densities and u ~ m v . This attraction is due to the N* resonance above the 
rjN threshold. Indeed, since the vacuum masses of rj, N and N* satisfy m v < 
mjy. — iriNi t ne rea l P ar t °f the self-energy has a negative value at uj ~ m n . At 
moderate densities when medium effects for the N and N* are not so strong 
and the inequality m rj < mjv*(p) — m^{p) still holds, the optical potential of 
the eta meson remains attractive for uj ~ m v . If some medium effects and/or 
a density increase bring the N* — h state down below the eta state, so that 
m v > win* (p) —m A r (p) j the rea l P ar t of the self-energy turns to be positive and 
the optical potential of the eta meson becomes repulsive pTUTll l|!12j . Therefore, 
the order of mesonic and particle-hole levels determines the sign of the self- 
energy, once two modes couple. Furthermore, the level difference between these 
modes is only about 50 MeV. Hence, the self-energy of the eta meson is very 
sensitive to the in-medium properties of the iV and iV* masses. In addition, the 
level coupling brings strong energy dependence to the self-energy. Therefore 
it may be wrong to discuss the eta meson in nuclear medium in terms of an 
energy-independent optical potential. In the next subsection, we discuss the 
spectrum of eta modes propagating in nuclear medium taking into account 
this energy dependence. 

The above discussion is quite general and does not depend on details of the 
mechanism of the N and N* mass reduction. 



2.2 Eta spectrum in the nuclear matter 

In this section we investigate the spectral function of the eta meson in the 
infinite nuclear medium. 
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Fig. 4. The same as in Fig. [3] but calculated with C = 0. 
The full in-medium propagator of the eta meson ([1]) can be decomposed as 



G v (u,k;p) = 



(7) 



where ct^'s are complex numbers being the solutions of the Dyson equation 
G~ x (o;, k; p) = with the self-energy ([3]): 
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- k 2 



mz 



9 2 P 



2u R (k]p) -iT N */2 



~aj-u R (k;p) +iT N */2] [u + u R {k; p)\ 



0. 



(8) 



We denote here u;i?(/c; p) = m* N ,(p)—m* N (p) + k 2 /(2m* N *(p)). This equation has 
two complex solutions with a positive real part, u>i(k) , i = 1, 2. These solutions 
describe two branches^ of the 77-meson spectrum in nuclear matteif^l. They 
are shown in Fig. [3] by solid lines for k = as functions of the baryon density. 
With the vanishing coupling constant g the eta-meson mode and iV*-nucleon- 
hole mode are decoupled and only the mesonic mode can carry the 77 quantum 
numbers. This case is depicted by the dashed lines on the left panel in Fig. [31 
The modes cross each other at p ~ 0.4 po for k = . The finite value of the 
coupling g leads to the avoided level crossing and the redistribution of the eta 
quantum numbers between both branches. 

The residua of the poles of the propagator G v (uj,k) at u = Ui quantify the 
occupation of each branch by the eta-meson quantum numbers: 



Ti(k;p) =2m ri Z i (k;p) , Zi(k; p) = ( 2u 



dll v (u,k; pY 
du 



■ (9) 



■Ji( k ) 



2 For the first time the two-branch spectrum of eta mesons in nuclear matter has 
been discussed in Ref. [56] . 

3 In general case the assignment of different solutions of the Dyson equation to the 
in-medium branches of particles and anti-particles should be done according to the 
sign of the real parts of the residua, cf. [57] 
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Fig. 5. The r\ meson spectral density (|10p in nuclear matter at densities O.lpoi 
0.5/3o an d lpo for the vanishing momentum. Dashed lines are calculated with the 
constant N* width Tn* = 75 MeV. The solid lines show the results with the energy 
and density dependent width shown in Fig. [2l The upper row is calculated with 
C = 0.2 and the lower one with C = . 

Note, it is important to include the imaginary part 7, by evaluating Z-factors. 
Then they have the following property ^sZ\{k) + QZ 2 (k) ~ and Ti(k; p) + 
T 2 {k; p) = 1 . The dependence of the occupation factors on the nucleon density 
is shown on the right panel in Fig.El We see that at p < 0.4 p the eta quantum 
numbers populate dominantly branch 2 and for p > 0.4 po branch 1. 

In order to illustrate the effects of the chiral symmetry restoration we depict in 
Fig.Hlthe same eta spectrum parameters as in Fig. [3]but calculated for C = 0. 
The mass gap of branch l,ui(k = 0), increases now with the increasing density 
and that of branch 2 decreases half as strong as for C = 0.2. The occupation 
factors show that branch 2 carries now the dominant part of the eta mesons. 

The spectral function is defined by the imaginary part of the full propagator 
(ID: 



cr(co>, k] p) = SsGJlj, k; p) . 

71 



(10) 



In Fig. the spectral density is shown for several values of the nucleon density 
at the vanishing momentum. At lower density (p = O.lpo) there is one peak, 
which corresponds to the eta mesonic mode. At higher densities a two-peak 
structure emerges. The second peak comes from the N*-h mode which couples 
to the eta mode in nuclear medium. The coupling strength increases as the 
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density increases. The dashed curves are calculated with the constant iV* 
width IV* = 75 MeV. The solid curves are calculated with the full energy 
and density dependent N* width depicted in Fig. [2J We note that dashed and 
solid curves coincide almost perfectly. For C = the low-energy peak of the 
spectral density is stronger pronounced than for C = 0.2 though it is located 
at somewhat higher energies. It is interesting to see how these differences in 
the spectral function reveal themselves in the finite nuclear systems and in eta 
production reactions. 



3 Eta meson modes in finite nuclei 

3.1 Density- averaged spectral function 

In eta production reactions on finite nuclei the spectral function ffTOl) is probed 
at different densities which produce superposed contributions to the observed 
final signal. Therefore, in order to get a feeling on the response of a nucleus 
to an external source with the eta-meson quantum numbers its instructive to 
study the spectral function averaged over the density distribution in a nucleus: 

a(u, k) = -jjdrr 2 p(r)Q G v (u, k; p(r)) . (11) 

Here p(r) is a spatial distribution of the nucleon density in a nucleus nor- 
malized to the total number of nucleons A = Att J drr 2 p(r). We use the 
2-parameter Fermi distribution: 

P(r) = - P \znS- (I 2 ) 

1 + exp (r_JS) 

For 11 B, which is the remaining nucleus for eta production reactions on 12 C, we 
have the central density po = 0.173 fm -3 , the nuclear radius R = 1.18A 1 / 3 — 
0.48 = 2.14 fm and the diffuseness parameter a = 0.5 fm. 

In Fig. [6] we plot the density-averaged nuclear response functions (JTTJ) with 
C = 0.0 and C = 0.2. Recalling the in-medium spectral function shown in 
Fig. G3> we find the following correspondence between the structures in these 
two response functions: the pronounced cusp structure at the eta production 
threshold corresponds to lower densities, when the eta is produced at the 
surface of the nucleus, and medium effects are small as seen on the left panel 
in Fig. [51 Two bumps — one in the continuum region (E > 0) and the other 
one in the bound regions (E < 0) — correspond to the two modes in the 
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Fig. 6. The density-averaged nuclear response function (jlip at zero momentum 
as a function of the eta meson energy, u, (lower scale) and the binding energy 
E = (uj 2 — m 2 )/2m v (upper scale) calculated for C = 0.2 (left panel) and C = 
(right panel). Dashed lines correspond to the constant N* width and solid lines to 
the energy and density dependent width. The vertical dotted lines mark the bare 
eta mass value. 



spectrum of the eta meson produced in the nucleus interior. The difference 
between the C = and C = 0.2 cases observed in Fig. [5] can be clearly 
seen also in the density-averaged response function. For C = 0.2 the bump in 
the continuum region is more pronounced and shifts to the higher energies, 
whereas the bump in the bound region slides down to lower energies and 
becomes smaller. Two modes repel each other stronger. Therefore, the response 
strength in the continuum region larger than that in the bound region can be 
a signal of the eta and iV* — h level crossing in the nucleus. 



It is worth noting that Eq. (TTTi) is not applicable for calculation of the bound 
state parameters, since the Green function calculated in the infinite nuclear 
matter satisfy the boundary conditions different from those for the Green 
function in a finite nucleus. Therefore, the infinite-matter Green function does 
not have correct analytic properties as the actual in-nucleus Green function 
and direct identification of the bumps in Fig. [6] with the poles of the Green 
function is not possible. We note that identification of the density-averaged 
spectral function ({TTI) with a local approximation of the actual in-nucleus 
Green function is not appropriate in regions where the nuclear density changes 
suddenly, like on the surface of a nucleus. Therefore, the cusp structures we 
observe in Fig. [6] can be only qualitatively compared to the actual spectra. 
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3.2 Bound and resonance states of eta-nucleus system 



In order to understand better the bump structure of the response function 
shown in Fig. 0, we have to study poles of the eta-meson-nucleus scattering 
amplitude on the complex energy plane. The positions of the poles correspond 
to energies and widths of bound states of the eta and N*-h modes in the 
nucleus. In practice, we search for the energy eigenvalues of the Klein-Gordon 
equation in the momentum space [60] for the eta meson with the angular 
momentum L in the nucleus: 

oo 

pH(p) + f p' 2 dp'v L (p,p';uj)^(p) = (oo 2 - mj) *(p) , (13) 



where \I/ is the wave function of the eta meson, and to is the energy of the eta 
meson. The non-relativistic binding energy is defined as E = (oo 2 —m 2 ) / (2 m v ). 
The potential vl(p,p'; oo) in the momentum space is given by the Fourier trans- 
form of the energy-dependent self-energy U v (oo,0] p(r)) evaluated in Eq.Q 
under the local density approximation: 

oo 

2 r 

v L (p,p';uj) = - dr r 2 11^,0; p(r))j L (pr) j L (p' r) (14) 

7T J 




with the spherical Bessel function jl(x). We put k = since we consider slow 
eta excitations created for the recoilless kinematic conditions. The boundary 
conditions can be implemented more easily by considering the Green functions 
in Lippmann-Schwinger equation rather than the wave function in the Klein- 
Gordon equation. The corresponding Lippmann-Schwinger equation written 
as an operator equation in the momentum space is 

* = ^(w)y(w)*, (is) 



where the Green function and the potential V are given in the operator 
form by G^(oj) = (uj 2 — m 2 — p 2 + it)" 1 and V(uS) = Vl(p,p';oj)p' 2 . The 
explicit form of Eq. (1151) is 



*(p) - / # : Lip f^,c *w = ° ■ <«) 

uj 2 — m l — p' 2 + it 



The bound state energies E, for which non-trivial solutions of the Lippmann- 
Schwinger equation (TT5l) exist, follow from the secular equation 



det (l - V(y/m* + 2m v E) G^(^m 2 v + 2m,£))=0. (17) 
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Some technical details of the calculations are given in Appendix |A] 

We discuss first the results obtained for L = and C = 0.2. On the first 
Riemann sheet we find two poles at E\ = (—88.53 — i 11.54) MeV and E 2 = 
( — 74.13 — i 15.00) MeV, which correspond to the deepest bound r\ state in the 
11 B nucleus and its first radial excitation, respectively. In the response function 
in Fig.[6]they manifest themselves as a small bump at about —80 MeV. On the 
first Riemann sheet there is another pole at E — (6.40 — i 0.46) MeV. Located 
above the threshold, this pole stands for a virtual state, which is not seen in the 
spectrum as a resonance bump but causes the cusp at the threshold. On the 
second Riemann sheet there is a resonance state at E — (61.02 — i 25.82) MeV, 
which shows up as a bump at ~ 50 MeV in the continuum region. Thus we 
conclude that the structure of the response function in Fig. [6] (left panel) 
contains signals of the two bound, one virtual and one resonance states. For 
the C = case the analytical structure changes. The attraction weakens 
and there remains only one bound state at (—22.18 — i 12.88) MeV on the 
first Riemann sheet. The other former bound state moves up and becomes a 
resonance at (3.47 — i 54.78) MeV on the second Riemann sheet. The other 
resonance state on the second Riemann sheet is located now at somewhat 
smaller energy (31.91 — i 54.79) MeV. The virtual state on the first Riemann 
sheet stays almost put at (6.13 — i 0.69) MeV. 



4 Formation spectra of rj mesic nuclei 

In this section we discuss the eta meson spectra in finite nuclei which can be 
observed in actual experiments. We consider nucleon pickup processes of the 
eta production on a nucleus. The (vr,p), (d, 3 He) and (j,p) processes on nuclear 
targets have been already proposed as formation reactions of the eta-nucleus 
systems. The advantage of the nucleon pickup processes is that it is possible 
to create the eta meson in nucleus almost at rest (recoilless production). The 
momentum transfer can be reduced to zero by adjusting the incident particle 
energy. For the recoilless production, we have a selection rule for the total 
angular momentum of the produced meson and the nucleon-hole state inside 
the nucleus. Since at the recoilless kinematics no angular momentum is brought 
to the residual nuclear system, the substitutional states are largely enhanced. 
As the result the contributions of the different angular momenta among the 
mesonic and nuclear states are strongly suppressed. For the 12 C target case, the 
dominant contributions are provided by the (0si/ 2 )~ 1 <S> s v and (0p 3 / 2 ) p 1 ®Vt) 
configurations, in which (nLs)~ l <8> £ v denotes ?7-hole configuration with the rj 
meson in the orbit I and the proton hole in the orbit L with the spin S and 
the principal quantum number n, although the ground state configuration is 
(0^3/2)^ 1 ® s v , which has different angular momenta for the rj and proton-hole 
states [7J. 
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In the experiments one observes the ejected nucleon in the final states inclu- 
sively and obtains the double differential cross sections d 2 <r/(df2d-E7v) as a 
functions of the energy, E^, and the emission angle of the outgoing nucleon. 
Investigating the shape of the missing mass spectra, we may hope to extract 
the in-medium properties of the eta and N*-h modes. Observing additional 
particles in the final state would help to confirm the creation of the eta and 
N*-h states in the nucleus (see, e.g., proposal [16J). 

Let us consider the 12 C(7,p) 11 B reaction as an example. The (d, 3 He) reactions 
were discussed in Refs. [TflOfllj and the (n,p) will be investigated in Ref. [58] . 
In the distorted-wave impulse approximation (DWIA) for the elementary pro- 
cess 7p — > rjp, the differential cross section is given by the product of the 
elementary cross section and the nuclear response function S(E): 



, > \ , , \ lab 




x S(£ 7 - E 



N 



where _E 7 is the incident-photon energy. In ([18]) we assume that the two-body 
process 7p — > rjp dominates the eta production and the direct production 
of the N*-h mode by the incident photon via the reaction jNN —>■ NN* is 
sub dominant. 



We calculate the nuclear response function S(E) following the Green function 
method for the hadrons in nuclei formulated in Ref. [55]. There the response 
function is given by the transition amplitudes T/(r) of 7 — > p~ l prj at the 
position r and the Green function G(E\ r, r') of the eta meson in the residual 
nucleus: 



S(E) = — Im J drdr'J2T}(r)G(E- AE jp ;r,r')T f (i 



(19) 



where the summation is taken over the all possible finial states and AEj is 
the hole separation energy measured from the ground state of the daughter 
nucleus. After the summation of the nucleon spins the transition amplitude 
Tf is given by 



T 



l/2,m s 



JM 



(20) 



with the proton-hole wave function ipj p , the distorted waves Xi an d Xj °f the 
incident photon and the ejected proton, the eta meson angular wave function 
Yi (r) and the spin wave function Ci/2,m a of the ejected proton. JM denotes 
the total angular momentum of the rj and proton-hole system. The Green 
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function G(E; r, r') is written formally as 



G(E;r,r' 



1 



$(r> _1 > 



(21) 



E - H„ + ie 



with the Hamiltonian containing the eta meson optical potential 



U v (uj,r) 



1 



U v (u; p(r)) . 



(22) 



2 m, 



In (l2~Tj) 0^ is the meson creation operator and is the residual nucleus 

state with a proton hole. The Green function G(E;r, r') given in Eq. (|2ip is 
calculated by solving the Klein-Gordon (or Schrodinger) equation with ap- 
propriate boundary conditions. In general, the Green function can depend on 
the nuclear states and give transitions between different states. Here we have 
assumed that these state dependence and transitions are negligibly small. 

The density averaged spectral function (TTTj) is related to the actual response 
function (fl9j) in the following way: since at the recoilless production the eta 
meson created in the nucleus moves slowly it cannot propagate a long dis- 
tance being absorbed rather soon. In this case the in-nucleus Green function 
G(E; r, r') given in Eq. fjlTLj) can be approximated as 



where G V (E; p(r)) is the Green function for infinite nuclear matter obtained 
in the local density approximation. Neglecting the distortion effects on the 
incident photon and the ejected proton in the transition amplitude Tf(r) in 
Eq. (l2"Ul) and replacing the product of the hole wave-functions, (r)ipj p (r) by 
the density distribution p(r) after the summation of the final states, we obtain 
an approximate response function 



This is the density weighted in-medium Green function introduced in f ill I) 
except for the normalization. According to the above argument, the approxi- 
mate spectral function is close to the in-nucleus spectral function obtained by 
solving the Klein-Gordon equation, if the local density approximation is good 
enough. In the proceeding calculations of the formation spectra of eta mesic 
nuclei we use the complete in-nucleus Green function f l2~Tj) . 

The calculated formation spectra of the eta meson in the (7,p) reaction with 
the 12 C target are shown in Fig. [7J The incident photon energy E 1 = 950 MeV 



G(E;r,r' 



)~5(v-v')G v (E;p(r)), 



(23) 




(24) 
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Fig. 7. Formation spectra of the eta meson in the (7,p) reaction on the 12 C target 
at = 950 MeV as a function of the energy of the proton, E ex , emitted at zero 
angle. The proton energy is counted from the eta production threshold Eq. The left 
and right panels are the formation spectra calculated in the chiral doublet model 
(the naive model) with C = 0.2 and C = 0.0, respectively. 

is chosen so that the recoilless condition is satisfied for the bare eta mass. We 
plot the missing mass spectra in the energy range wider than those shown 
in Refs. [101)111)12] to see the global structure of the spectra and to be able 
to identify the bound and resonance states discussed in the previous section. 
The energy of the emitted proton is counted from the creation threshold of 
the vacuum eta mass. The left (right) side of the threshold corresponds to the 
bound (continuum) energy region. We have calculated the sum of all contribu- 
tions from the possible configurations of the proton-hole states, S1/2 and P3/2, 
and the eta states with I = to 6. The total spectra are dominated by two 
contributions from the (0si/ 2 )p 1 ® s v and (0p 3 / 2 )p 1 ®p v configurations due to 
the selection rule at the recoilless kinematics. The separation energies of the 
s state protons are taken into account in the spectra. The formation spectra 
are discussed in more details in Refs. [101)111)12] . 

The missing mass spectrum calculated in the chiral doublet model with C = 
0.2 and shown on the left panel in Figl7]has the richer structure than the spec- 
trum obtained with C = 0.0, when no level crossing takes place. We compare 
the sub-component (Osi^)^ 1 <8> s v in FigJ7|with the response function plotted 
in Fig. [6j The centrifugal potentials are not taken into account in the latter 
function. We find that still the approximated response function reproduces 
well the bump structure of the spectrum obtained in the full calculation with 
the actual Green function. (Exception is the eta threshold energies, when the 
contributions from the eta meson on the surface of the residual nucleus are 
dominated.) This is because the distortion factors have smooth energy depen- 
dence and the initial photon is less distorted. Thus, the conclusions drawn 
in the previous section for the approximated spectra hold also for the actual 
formation spectra. Namely, the bump structures correspond to the bound and 
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resonance states. For the C = 0.2 case, we have found the two bound states 
and one resonance state. They correspond to the small bump around E = —80 
MeV and the pronounced peak in E = 80 MeV, respectively. 

We conclude that the more pronounced bump structure in the continuum 
energy region can be a signal of the level crossing. Note that the recoilless 
kinematics largely suppresses the second bound state (the radial excited state) 
in the spectrum, since there are no proton-hole states in ll B to satisfy the 
selection rule with the second s state of the eta meson. 



5 Conclusion 

We have investigated the eta-nucleus system, paying attention to the iV(1535)- 
nucleon-hole mode in nuclear medium. We have found that the N*-h contri- 
bution is important for the eta meson in nuclear matter and induces strong 
energy dependence of the eta self-energy. We have also shown that reduction 
of the difference between the masses of N* and iV as a result of the medium 
effects leads to level crossing of the eta and N*-h modes, and modification of 
the eta spectral function. The change of the spectral function is clearly seen 
in formation spectra of eta nuclei. 
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A Appendix 

The numerical solution of the equation (fl6|) is performed by the method sug- 
gested in Ref. [50] • The integral equation is written as a matrix equation in 
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momentum space by replacing momentum integral with a sum over N grid 
points: 



N 



ty m = J2M mn V n (A.l) 



n=l 



with 

M mn = V mn (uj)G^(u)W n , (A.2) 
where W n is a weight function for the discretization of the integral and 



V mn (u) = v L (p m ,p n ;u)p n , (A. 3) 

1 

uj' a — vri l r] — p% + ie 

At the end we search for zeros of the determinant of the N x N matrix: 



Gi 0) H= , „, _ 9 , ... . (A.4) 



det(l — M) = (A.5) 



for the complex binding energy E = (u 2 — m 2 )/(2m v ) . 

In this method, one can choose properly the boundary conditions of the scat- 
tering equation. For the first Riemann sheet solutions, we just solve Eq. (1A.5l) 
for a complex E. For the second Riemann sheet, where we can obtain reso- 
nance solution, one should treat ie in the Green function given in Eq. (1A.4I) 
before going to the complex plane. This can be done as follows. The integral in 
the left hand side of Eq. f|T6l) can be written in terms of principle value integral 
as 



*(p) = P I dp' ^^^ ip') ~ 4 *kv L (p, k; coMk) , (A.6) 
J uj l — — p' 1 2 



where k is a complex momentum k 2 = u 2 — m 2 = 2 m v E on the second 
Riemann sheet. As discussed in Ref. (61], these terms are implemented to the 



matrix formulation by extending the matrix M in flA.5j) and defining an extra 
grid point Pn+i = k as 



M N+1>n = v L (k,p n ,uj) , ( A.7) 



WnPl 
k 2 — p 

/ N -W m k 2 1 \ 
M U)N+ i = v L (p n , k; u) I 53 u2 m 2 ~ i o 7rk ) > ( A - 8 



^k 2 -p 2 m 2 
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/ N -W m k 2 1 \ 
M N+ljN+1 = v L (k, k;uj)[J2 u2 ~ 2 _ ■ ( A - 9 ) 

\ m =l K -Pm A ) 
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